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I. INTRODUCTION 



Three-nucleon forces (3NFs) are an indispensable ingredient in accurate few-nucleon and nuclear structure calculations. 
In particular, the three-nucleon system shows clear evidence for 3NFs, see the recent general introduction [l| . Chiral 
' Cji ' effective field theory is the appropriate tool to analyze nuclear forces. Precise two-nucleon potentials have been 
^ \ developed at next-to-next-to-next-to-leadingorder (N'^LO) in the chiral expansion, see Refs.JpJS- 3NFs first appear 
Ch . at N^LO in the Weinberg counting scheme P, [E] and have been analyzed and scrutinized in [a ItT sl lol. [Tol ITU . There 
are various reasons to improve the theoretical precision of these 3NFs: 1) one should utilize the two- and three-nucleon 
, forces at the same order in the expansion, 2) there are some outstanding discrepancies between theory and experiment 
^ ' like e.g the recently measured differential cross section in deuteron break-up at low energies or the long-standing 
! Ay puzzle [T^ . [T^ and 3) the theoretical uncertainty employing only the leading 3NFs quickly grows with increasing 
energy if one investigates e.g. nucleon-deutcron scattering or break-up. It is therefore timely and necessary to derive 

■ the chiral 3NF at N^LO. In this work, we focus our attention on the long-range contributions at this order, which are 
. • free of unknown coupling constants. Here, we consider the chiral effective Lagrangian with pions and nucleons. The 

precise relation of the results presented here to an effective field theory including also spin-3/2 degrees of freedom will 

■ be the subject of a subsequent paper. 

^p. I As we will show later, there are five different topologies contributing to the 3NF at N^LO. From these, three topologies 
^ . make up the long-range contribution, which is defined by not including any multi-nucleon contact interactions. This 
long-range part is given by two-pion exchange (27r) graphs, two-pion-one-pion exchange (27r-l7r) graphs and the so- 

■ called ring diagrams, where the pion loop connects all three nucleon lines. The first type of graphs has recently 
d ' been considered based on the so-called infrared regularization in Ref. p^ . we will compare our results to that work 

below. The spin-isospin structures originating from ring diagrams with one explicit delta intermediate state have 
been incorporated in the Illinois 3NF model [16]. Earlier, Coon and Friar (l7| systematically constructed the 1/m 
(with m the nucleon mass) corrections to the 27r exchange 3NF, and the so-called drift effects due to the boost of the 
two-nucleon force were worked out in (isf . We will compare to these works in our subsequent paper, where we discuss 
the shorter-ranged part of the 3NF at this order and the corresponding 1/m corrections. 

This manuscript is organized as follows. In Sect. |TT]we first write down the effective chiral Lagrangian and discuss 
the general structure of the 3NF at N'^LO. In the subsequent subsections III Al III Bl and III Cl we give our results for 
the 2tt, 27r-l7r and ring graphs, respectively, both in momentum and coordinate space representations. We end with 
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a summary and conclusions in Sect. IIIII In particular, we make some comments on the 3NF arising in an extension 
of this work with explicit delta degrees of freedom. The lengthy expressions for the momentum space representation 
of the ring diagrams are relegated to the appendix. 
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II. LONG-RANGE CONTRIBUTIONS TO THE 3NF AT N^LO 



The calculations performed in the following are based on the effective chiral Lagrangian for pions and nucleons. We 
employ here the heavy baryon formulation and display only the terms of relevance for our study: 

= £^=" + £^=1 + £^=2 ^ ^ ^ ^ 

^ ^{y^UV^U^ +X+)+N{^vD + gAU-S)N + ... , 

= N (ci (x+) + C2 (w • u)' + C3 u • w + C4 [S^, S'']Uf,u,) N + ... , 

= N (di65 • u{x+) + idisS^'iD^, X-] + d28{^{x+)v ■ D + h.c.) + . . .) TV + . . . , (2.1) 

where the Ci are low-energy constants and N , and 5*^ denote the large component of the nucleon field, the nucleons 
four- velocity and the covariant spin vector, respectively. We use standard notation: U{x) = v?(x) collects the pion 
fields, Uf^i = i{u^d^u — ud^u''), x+ — u'^xu'' + ux''u includes the explicit chiral symmetry breaking due to the finite 
light quark masses, (. . .) denotes a trace in flavor space and I?^ is the chiral covariant derivative for the nucleon field. 
Notice further that the first terms in the expansion of [/(tt) in powers of the pion fields read: 

, i 1 9 , , Q 8q! — 1 i , , 

t^w = i + :^^-^-^^ -;^(^-^) +^:^^ +••• ' (2-2) 

where r denote the Pauli isospin matrices and a is an arbitrary constant. For further notation and discussion, we 
refer to Ref. [l^ , a recent review is given in Ref. [l^l . Following Weinberg [1, Q , the dimension A of the Lagrangian 
is defined via 

A = fi+^7i-2, (2.3) 

where d and n are the number of derivatives or insertions of the pion mass Mj^ and nucleon field operators, respectively. 
The pertinent low-energy constants (LECs) of the leading-order effective Lagrangian are the nucleon axial-vector 
coupling gA and the pion decay constant F^^. Notice that while all couplings and masses appearing in the effective 
Lagrangian should, strictly speaking, be taken at their SU(2) chiral limit values, to the accuracy we are working, we 
can use their pertinent physical values. In addition, we have the LECs dig, dig and ^28 from the TriV Lagrangian at 
order A = 2. The ellipses in the brackets in the second line of Eq. (j2.H) refer to terms proportional to the LECs 
di, 2, 3, 5, 14, 15 and ^24, 26, 27, 28, 30 which generate mrNN vertices [19J but do not contribute to the 3NF at N'^LO as will be 
shown later. We also omit in Eq. (|2.ip pion vertices with A = 2 and proportional to LECs which do not show up 
explicitly in the formulation based on renormalized pion fields at the considered order, see Ref. [2l| for more details. 

For a connected A'^-nucleon diagram with L loops and Vi vertices of dimension Ai, the irreducible contribution^ to 
the scattering amplitude scales as Q'^ , where Q is a generic low-momentum scale associated with external nucleon 
three-momenta or and 

= -A + 2N + 2L + Y^ VAi ■ (2.4) 

i 

Consequently, at N'^LO, which corresponds to = 4, one needs to take into account two classes of connected diagrams: 
tree diagrams with one insertion of the A = 2-interactions and one-loop graphs involving only lowest-order vertices 
with A = 0. Notice that it is not possible to draw 3N diagrams with two insertions of A — 1-vertices at this order. We 



^ This is the contribution which is not generated through iterations in the dynamical equation and which gives rise to the nuclear force. 
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(a) (b) (c) (d) (e) 

FIG. 1: Various topologies that appear in the 3NF at N^LO. Solid and dashed lines represent nucleons and pions, respectively. 
Shaded blobs are the corresponding amplitudes. The long-range part of the 3NF considered in this paper consists of (a) 2ii 
exchange graphs, (b) 27r-l7r diagrams and the so-called ring diagrams (c). The topologies (d) and (e) involve four-nucleon 
contact operators and are considered of shorter range. 



further emphasize that similar to the case of the leading four-nucleon force considered in Refs. [22, llj], disconnected 
diagrams lead to vanishing contributions to the 3NF and will not be discussed in what follows. 

The structure of the 3NF at N'^LO is visualized in Fig. [T]and can be written as 

— + ^2v-lTT + ''ring + '^iTr-cont + ''27r-cont + ^l/m ' \'^-^) 

While the 27r-l7r, ring and two-pion-exchange-contact (27r-cont) topologies start to contribute at N^LO, the 27r and 
onc-pion-exchange-contact (l7r-cont) graphs already appear at N^LO yielding the following contributions to the 3NF 



(3) _ 9a CTl ■ gl g3 ■ g3 



Ti ■ T3 (-4ciM^ + 2c3 qi ■ 93) + C4T1 X T3 ■ r2 gi X 53 • ctj 



x.(3) 9aD as-qs ^ ^ 

* ^ -jFf ■ "3 "1 ■ ' (2.6) 

where the subscripts refer to the nucleon labels and qi = pi' — pi, with pi ' and pi being the final and initial momenta 
of the nucleon i. Further, qi = |qi|, ai denote the Pauli spin matrices and D refers to the low-energy constant 
accompanying the leading ttNNNN vertex. Here and throughout this work, the results are always given for a 
particular choice of nucleon labels. The full expression for the 3NF results by taking into account all possible 
permutations of the nucleons^, i.e.: 

Vgf^ii = V3N + aU permutations . (2.7) 

In this work, we focus on the long-range contributions V2^\ ^2T-i7r ^^'^ ^ing resulting from diagrams (a), (b) and 
(c) in Fig. [Tl respectively. The remaining terms resulting from graphs (d) and (e) in Fig. [1] and the relativistic 
1/m-corrections will be discussed separately. 



A. Two-pion exchange topology 



The 2tt exchange diagrams at v = 4 corresponding to topology (a) of Fig. [T] are depicted in Fig. [51 We do not show 
graphs involving tadpoles at pion lines which just renormalize the pion field and mass, see [2l| for details. A close 
inspection of diagrams in Fig. [2| shows that most of them do not generate 3NFs. First, the contribution from graphs 
(12-17) involves an odd power of the loop momentum I to be integrated over and thus vanishes. Second, diagrams 



^ For three nucleons there are altogether 6 permutations. 
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FIG. 2; Two-pion exchange 3N diagrams at N'^LO. Solid dots (filled rectangles) denote vertices of dimension Ai — [Ai — 2). 
Diagrams which result from the interchange of the nucleon lines and/or application of the time reversal operation are not 
shown. For remaining notation see Fig. [T] 



(1), (2), (4), (8), (10), (30-32), (34) and (35) just renormalize the external nucleon legs. Similarly, Feynman diagrams 
(3), (9), (22-24), (26-28) lead to renormalization of the leading pion- nucleon coupling without producing any new 
structures. All these contributions are taken into account by replacing the bare LECs in the leading 27r exchange 3N 
scattering amplitude by renormalized ones. This suggests that there are no N'^LO corrections to the 3NF from these 
graphs since the 27r exchange 3N diagrams at order = 2 do not generate any nonvanishing 3NF. Given the fact 
that nuclear potentials are, in general, not uniquely defined, the above argument based on the (on-shell) scattering 
amplitude should be taken with care. We have, however, verified that this is indeed the case by explicitly calculating 
the corresponding 3NF using the method of unitary transformation along the lines of Ref. [23| . From the remaining 
graphs in Fig. [21 diagram (11) does not contribute at the considered order due to the 1/m-suppression caused by the 
time derivative entering the Weinberg- Tomozawa vertex.'^ For the same reason, diagram (25) also leads to a vanishing 
result at the order considered. Here, the time derivative acts either on the pions exchanged between two nucleons 
leading to a 1/m-suppression or on the pion in the tadpole giving an odd power of the loop momentum to be 



^ This graph does not involve reducible time-ordered topologies. Its contribution to the nuclear force is, therefore, most easily obtained 
using the Feynman graph technique. The 1/m-suppression from the time derivative entering the Weinberg- Tomozawa vertex follows 
then simply from the four-momentum conservation. 
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integrated over. Further, it is easy to see that Feynman diagrams (18) and (21) do not contribute as well. Diagram 
(29) involves one insertion of the ttttNN vertices of dimension v — 2. The relevant vertices are proportional to LECs 
'^1,2,3,5,14,15 ^-iid (^24, 26, 27, 28, 30- The Corresponding 3NF is shifted to higher orders since all these vertices involve at 
least one time derivative, see [l9l | for explicit expressions. Last but not least, we also found that diagram (33) does not 
generate any 3NF. Thus, we are left with diagrams (5-7), (19) and (20). The 3NF contribution from diagrams (5-7) 
can be evaluated straightforwardly using the expressions for the effective Hamilton operator from Ref. [2^ . Diagrams 
(19) and (20) do not involve reducible topologies and can be evaluated using the Feynman graph technique. Notice 
that the individual contributions from graphs (19) and (20) in Fig. [5] and from diagram (20) in Fig. [3] depend on 
the arbitrary constant a which specifies the parametrization of the matrix U, see Eq. (|2.2p . Clearly, their sum is 
a- independent. 

We are now in the position to present our results. The expressions for diagrams (5-7) and (19) can be cast into the 
form of Eq. (|2.6p leading only to shifts in the values of the LECs c; 

9a - 9a - 9\M^ 

with 6ci = -0.13GeV"i and 5c^ ^ -(5c4 = 0.52 GeV-^ These shifts are of the order of 20% to 30% of the 
corresponding LECs, and thus can not be neglected in precision studies of 3NFs. In contrast to this, the contribution 
from graph (20) takes a more complicated form compared to Eq. (|2.6p and is given by 



XI • X3 {M^{Ml + iql + iql + 4<fi • 93) -f {2MI + ql+ql + 2qi ■ qs) (2.9) 



- 256^F6 [ql + A'miqi + M^] 

X (3M2 + 3ql + 3ql + Aqi ■ q3)A{q2)) - Ti X T3 ■ qi X q^ ■ (M, + (4Af2 + 9? + g| + '2q^ ■ 93)^(92)) 

Here, we have used dimensional regularization to evaluate the loop integrals. In this framework, the loop function 
A{q) is given by: 

A(9) = ^arctan^. (2.10) 

We further emphasize that the above expressions correspond to the choice a — 0. Notice also that some pieces in 
Eq. (12. 9p can be brought into a form corresponding to the 27r-l7r and 27r-cont topologies by canceling out the pion 
propagators with terms in the numerator. The 27r exchange contributions arising from diagrams in Fig. [5] have also 
been considered recently by Ishikawa and Robilotta based on the infrared regularization 15]. We have verified that 
the long-range 2tt exchange contributions in Eqs. (|2.6p . (|2.8p and (|2.9p agree with the ones given in Ref. [3l provided 
the latter are expanded in powers of 1/m and the leading terms are kept. 

The coordinate space representation of the 2tt exchange 3NF can be obtained in a straightforward way. For the leading 
terms in the first line of Eq. (|2.6p and the corrections in Eq. (|2.8p one gets 

V.An.,r,,) = /|^|^e--e--y..(9l,93) 



5^m;; ^ - ^ - 

(71 • V12 (73 • V32 



Ti • T3 4ci + 2C3V12 • V32 + C4T1 XT3 ■T2 V12 X V32 • CT2 



1287r2F4 

X C/i(xi2) Ui{x32), (2.11) 

where fij = fi — fj is the distance between the nucleons i and j, the Xi = M^^ ri are dimensionless distances, the 
act on Xi and Xij = \xij\. Further, the scalar function Ui{x) is defined as: 

Similarly, one obtains for the terms in Eq. (|2.9p which do not involve the loop function A{q2): 



V2^(ri2,r32) = 77T?d^^ • V12 (73 • V32 X1-T3 (-l + 3V22+3V22+4Vi2-V3 



9\Ml 
40967r3F6 

Ti X r3 • r2 V12 X V32 • CT2 Ui{xi2) Ui{x32) ■ (2.13) 
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FIG. 3: 27r-l7r diagrams at N^LO. Graphs resulting from the interchange of the nucleon Unes are not shown. For notation see 
Figs. [U [2 



It should be understood that the obtained expressions are only valid in the region of space where the interparticle 
distances are large (i.e. larger than the inverse pion mass). The behavior of the potential at shorter distances is, in 
general, affected by the regularization procedure which is not considered in the present work. 



To obtain the coordinate space representation for the terms in Eq. (|2.9p which involve the loop function A{q2), it is 
more convenient to proceed in a different way in order to avoid a complicated angular integration: 



1^27r(ri2, ^32 ) 



d^qi (fiq2 d^qs ,^,3.3 



(27r)3 (27r)3 (27r)3 



i2nyS'iq'i + 92 + 93 ) e*«^"-^ e^^^"-^ Vs^gl, 9*2, 93) 



d'ro 



d^qi d^q2 d^q3 



^iqi-rio „iq2-r20 piqz-rzo 



(27r)3 (27r)3 (27r)3 



V27r(<?l, (72, 93) 



■ CTi • V12 (?3 • V3 



Ti -xs (2- 



40967r3F6 

+ T-l X T3 • T2 V12 X V32 • CT2 (4 - V?2 - V32 - 2Vi2 • V32) 

X ^ / d'^xUi{\xx2+x\) Wi{x) Ux{\x32+X\), 



2 - V^2 - 2Vi2 • V32)(3 - 3V^2 - 3V^2 - 4Vi2 • V32) 



where 



Wi{x) 



d3 



M2 J (27r)3 2x2 ■ 

For various techniques to evaluate the integral in the last line of Eq. (|2.14p the reader is referred to Ref. [if 



(2.14) 



(2.15) 
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B. Two-pion— one-pion exchange topology 



Consider now the 27r-l7r 3NF arising from diagrams shown in Fig. [3] corresponding to topology (b) in Fig. [T] They 
can be written in the form: 



V2tt-1tt 



^3 ■ g3 

+ CTi • 93 Fr^iqi) + iT2 ■ gi Fe{qi) + 0^2 • 93 Fjiqi)) + ti x T2 • T3 cti x (T2 • qi Fs{qi) 



Ti ■ Ta (ct2 ■ qi qi ■ 93 ^1(91) +^2- qi F2{qi) + ^2 ■ qs. F^iqi)) + T2 ■ T3 (cti • <fi (fi • 93 -^4(91) 

(2.16) 



with Fi....,8('?i) being scalar functions. Using the formal expressions for the effective Hamiltonian from Ref. [23| we 
obtain the following result for the first six graphs: 



Fi{qi) 
Fsiqi) 
Fi{qi) 



9a 



2567rF6 



2M, 8M2 + ql 



4Afi 



ql 



ql 



ql 



A{qi 



9a 



2567rF6 
4^5(91) 

qi 



3Af, + (8A/2 + 3q2)A(gi) 
A(<zi). 



9\ 



1287rF6 



(2.17) 



Consider now diagrams (7-15) which involve one insertion of the Weinberg- Tomozawa vertex. The first three graphs 
only contribute to the functions -^1.3(91): 



^^i(<Zi) 
F:i{qi) 



9\ 



2567ri^6 



9\ 



. 9i 

M. 



ql 



{ql + AMl)A[q^) 



(2.18) 



2567ri^6 

We find that diagrams (10) and (11) do not generate a 3NF while the contribution from graphs (12-15) has the form: 



F2{qi) 



9\ 



12877^6 



+ {ql + 2Ml) A{qi 



(2.19) 



The next two diagrams (16) and (17) do not involve reducible topologies and can be dealt with using the Feynman 
graph technique. It is easy to see that their contribution to the scattering amplitude is suppressed by 1/m due to the 
time derivative which enters the Weinberg- Tomozawa vertices. The next two diagrams (18) and (19) yield vanishing 
contributions to the 3NF for exactly the same reason as do the two-pion exchange diagrams oc g'^ at order v = 2. 
Next, the contribution from diagram (20) reads: 



F6(gi) = 2^^7(91 
Fs{qi) - 



9\ 



9\ 



647rF6 



A'U + (ql + 2A'C)A{qi: 
{ql+4M^)A{q,) . 



(2.20) 



5127ri^6 

Finally, it is easy to see that the last Feynman diagram (21) in Fig.[3]yields a vanishing result at the order considered. 

The coordinate space representation of the 27r-l7r 3NF can be obtained straightforwardly employing the same Fourier- 
type integrations as in the the first line of Eq. (I2.14p . This leads to the following expression where short-range terms 
resulting from constant contributions to ^2, 3, 6, 7. 81 which are proportional to M^r, are not shown: : 



^27r-l7r(?'l2, ^52 ) 



9\Ml 
8192 TT^F2 



0^3 • V32 2ti • T3 



(72 • V12 V12 • V3 



2glUi{2xi2) + [l + 9l)Wi{xi2) 



+ (1 - 2gl)W3{xi2)) - 2(72 ■ Vi2(2M^i(a;i2) -I- 1^2(2:12)) - <?2 • V32(4(l - 2g\)Wi{xi2) 



+ {1 ~ 3gl)W2{xu)) 



- 4t2 • T3 



9a ^1 • '^12 V12 • V32 Wi{xi2) +g\ai- V32 ^^2(2:12) 



+ 2 CT2 ■ V12 {2Wi{xi2) + 1^2(2:12)) + CT2 ■ V32 (2M^i(a;i2) -I- ^^2(2:12)) 

n X X2 • T3 (71 X CT2 • V12 (4M^i(xi2) + ^^2(2:12)) ) Ui(x32) , 



(2.21) 
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FIG. 4: Ring diagrams at N'^LO. Graphs resulting from the interchange of the nucleon lines are not shown. For notation see 
Figs. [11 [2] 



with 



and 



W2{x) = -ViWi{x)^—^{l + 2x{l + x)), 

x^ 



W^ix) 



M^J (27r)3 



A{q) 



2Ei{-2x) 



-2x 



(2.22) 



Ei{x) = - 



t 



(2.23) 



We emphasize again that the above expressions are only vahd at large distances. 



C. Ring diagrams 



We now regard ring diagrams shown in Fig. |4] which correspond to topology (c) in Fig.[TJ These are most cumbersome 
to evaluate. The contributions from the first two diagrams can be obtained using the expressions for the effective 
Hamilton operator given in Ref. [1^. This leads to the following structures: 



* ring 



+ 



* ring 



UJa LUb 

4 



U)^ UJb UJc UJa UJc 



(2.24) 



UJa UJb UJ-^ 

where represents the spin, isospin and momentum structure which results from the vertices entering the diagram 
i and uj denotes the pion free energy, uj = \J + . Substituting the expressions for the vertices, the result can be 
written in the form 



with the numerator 



1 



(27r)3 



cPh cPh cPh S'^ih -h- qi) S'^ih - h 



qs) 



API 



(2.25) 



V — — 8ti ■ T2 h X I3 ■ (72 h X I2 ■ (?3 h ■ h 
+ 6I2 X Is ■ ai h X I2 ■ CT3 h ■ h ■ 



4ti ■ T3 h ■ I2 h -hh-h + 2ti x T2 ■ h x h ■ S2 h ■ h h ■ k 

(2.26) 



Carrying out the two trivial integrations over, say, li and I2 leads to the standard three-point function integrals. The 
latter can be evaluated but the resulting expressions are rather involved, see appendix |^ It is more convenient to 
evaluate Eq. (|2.25p in configuration space using again the same definition as in the the first line of Eq. (|2.1ip . This 
leads to the following compact result: 



V;ing(rr2, ^52)= ( ^ 



(fh d^h d^h 



(27r)3 (27r)3 (27r)3 



[q + Am [q + MiYn + Mi] 



4096 7r3i^6 



4X1 • T2 V23 X V12 • (72 V23 X V31 • 1T3 V31 • V12 



2X1 ■ T3 V23 • V31 V23 ■ V12 V31 • V12 + Xi X X2 • X3 V23 X V12 • CT2 V23 • V31 V31 • V12 



3V31 X V12 • CTl V23 X V31 • CT3 V23 ■ Vis 



C^i(a;23) C^2(a;3i) Ui(xi2) , 



(2.27) 



where the derivatives should be evaluated as if the variables ^12, a?23 and X31 were independent^ and we have introduced 



U2{x) = SttM^ 



(2^ [(72 + M2]2 



(2.28) 



Consider now diagrams (3-5) in Fig. [3] which involve one insertion of the Weinberg- Tomozawa vertex. The correspond- 
ing contribution to the 3NF can be evaluated in a similar way as described above and using the explicit expressions 
for the effective Hamilton operator from Ref. |23)]. We find that graphs (3) and (4) lead to vanishing 3NF while the 
contribution from diagram (5) reads: 



K-ing(7^12, ^52 ) = 



9\Ml 



2X1 • T2 (V23 • V31 V31 • V12 - V31 X V12 • (Ti V23 X V31 • (73 ) 

Ul{x21,) Ui{x:ii) Ui{xi2) ■ 



2048 7r3 F6 

+ Xl X X2 • X3 V31 X V12 ■ (71 V23 ■ V31 



(2.29) 



Finally, we found that the contributions from the last three diagrams in Fig. [3] are suppressed by the factor 1 /m and, 
therefore, need not to be taken into account at the order considered. 



III. SUMMARY AND CONCLUSIONS 



In this work, we have derived the long-range contributions of the 3NF in chiral effective field theory at N^LO in the 
chiral expansion. These results constitute the first systematic corrections to the leading 27r exchange that appears 
at N^LO in the Weinberg counting. We have shown that there exists three different topologies that contribute to 
the long-range components of the chiral 3NF. While some of these simply renormalize the LECs of the leading order 
27r exchange topology, there are many new spin-isospin structures generated at N^LO from all three topologies. We 
have given the momentum and coordinate space representations of the various contributions. At this order, there is 
an arbitrariness in choosing the multi-pion interactions, cf. Eq. (|2.2[) . which allows to shuffle strength between the 
various contributions. The sum of all diagrams contributing to the 3NF is, as demanded by field theory, independent 
of the choice of this parameter. To be precise, we have given our results for the special choice a = 0. It is important 
to stress that the corrections to the 3NF discussed here are free of any unknown low-energy constants and can be 
expressed entirely in terms of the nucleon axial- vector coupling gA, the pion decay constant -FV and the pion mass 
M^. It should further be emphasized that the dimension- two LECs C3 and C4 which enter the leading 3NF and would 
also contribute at N^LO are larger than their natural values Ci ~ f;yi/(47rF^) ^ 1 GeV~^ [l^]. This is understood 
in terms of A(1232) s-channel and vector-meson t-channel excitations [2J]. It is, therefore, expected that the theory 
with explicit spin-3/2 degrees of freedom leads to an improved convergence as compared to the pure pion- nucleon 
theory considered here. Of course, the theory with explicit deltas leads to many more structures and has also been 
much less developed phenomenologically. Still, it is important to confront the results obtained here with similarly 
accurate calculations in the delta-full theory. Work along these lines is in progress. 
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APPENDIX A: EXPRESSIONS FOR RING DIAGRAMS IN MOMENTUM-SPACE 

In this appendix we give lengthy expressions for ring diagrams in Fig. [4] in momentum space. The contributions from 
diagrams (1) and (2) can be expressed as: 

Vring = CTl • (72 T2 • Xs Ri + Si ■ ^1(72 ■ qi T2 ■ T3 i?2 + CTl • gil72 ■ T2 ■ T3 i?3 + CTl • q3^2 ' <Zl T2 • T3 R4 

+ (?i • q3^2 ■ qs T2 ■ T3 i?5 + Ti • T3 Re + ai ■ qias ■ qi Rj + ai ■ qiaz ■ qs Rs + ■ q03 ■ <Zi -Kg (A.l) 

+ CTl • CT3 i?10 + gi • X (72 Tl • X Ts Rll. 

where the functions Ri = Ri{qi, qs, z) with z — qi ■ q^ are defined as follows: 

-1 + z2) g^M^ {2Ml + ql) {q^qs + AM^ (zq, + q^)) A {q2) g\ql (?! + zq^) + zq^ {-qf + qj)) 



Ri = 



USF^TT (4 (-1 + z2) M2 - g2) (4M2g3 + ql) 128^% (-1 + z2) qiql 

A (ga) g\ {zql jzqi - ga) q^ + 2Af2 (_2 + z^) qf - {l + z^) giga - zql)) 

128F67r(-l + z2)q^q3 

A {qi)g\ (2M2g2 + ^3 (_^g3 ^ (2 - iz^) qfq^ - z (-2 + z^) q^qj + qj)) 

128F'^Tr + z^) ql 

(8 {-1 + z ) {2zqi + (1 + z ) ga) + g2g3 [z gi + 



i?2 = 



32F6 (-1 + z2) (4 (_1 + ^2) ^2 _ ^2) 

z (-1 + z2) giga - g2) + 2AC {z (-2 + z2) g? - (l + 2z2) g^ga + 3z (-2 + z^) 91932 + (-3 + 2z4) gf)) 

A (g2) gjgl (-2Afg ((1 + z2) q, + 2zga) + zga ((l + z2) qf - 2g32)) ^ 
128F67r(-l + z2)2g3g2 

A (ga) 5i (^2 (2zg2 + (l + 3z2) q^q^ + 2zql) + zga (-zg? - z2g2ga + zgigf + gf)) , 



64Ffi7r(-l + z2)%3g3 
{2Ml ((1 + z2) 92 + ^ (3 + q^q^ + (1 + ^2) ^2^ ^ 



^'^>A/^2^^ I -,2\ „2 I ^ z Q I „2 

128i^6^(-l + z2)2 92g2 



R3 



ga (- (z + z3) g? + (2 - Sz^ + z^) g^ga + z (l + z^) 9192 ^ ^ ^2^ ^3^^ _ 

/(4 : 0,~gi,g3;0).g^ .4 / „ 2 2 , , 2\ 2\ 

32FeM + z2fg2(-4(-l + z2)M2 + ,2)^3 (g2g3 (-2z g^ + (l + z ) ga) - 

8(-l + z)(l + z)Af^ (z (2 + z^) g? + (1 + 2z2)%2q,^ + z (2 + 7z'^) gig| + (l + 2z2) qfj + 

2Mlql {2zql + (1 - z^ + 6z'^) qjqs - 2z (-1 - 3z2 + z'^) gig| + (3 + 3z^ - 4z'^) gf)) 

g^M^ (2M2 + 932) (g2g3 + 4jv/2 (^g^ + g^)) 
+ 128i^%g2 (4 (-1 + z2) A/2 - g2) (^^^^q^ + ^3) ' 

zA (ga) ff^gj (-4A/2 (g^ + zqs) + ga (2zg2 + (-1 + z2) q^q^ - 2zg§)) _ 
128F^7r {-1 + z^ f qfql 

^"'^ (M'^ {-2z (-3 + z^) g? + 4 (1 + z^) g^ga + 4zga^) + ga (- (l + z^) qf- 



128F^TT i-1 + z^f qlql 
2z3g?ga + (l + z^) qig2 ^ 2zg|)) - 

zA (gi) g6 (2M2 (2g2 + 4zgig3 + (l + z2) qj) + q^ {-2zqf + (l - 3z2) ^2^3 + 2zgig2 + (l + z2) ^33)) 



128i^%(-l + z2) Qigf 

/(4 : 0, -gi, gs; 0)z5^ / 4 //, , 2\ 2 , / , , 2\ , 2\ 2\ 

(g2g3 ((1 + 2 ) gi + z (-1 + z ) gigs - (l + z ) ga) 



32F6 (-1 + z2)" gi (-4 (-1 + z2) Af2 + g2) ^2 
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8(-l + z)(l + z)Af4 {^Szcil + (-1 + lOz^) qlq^ + 3z (l + 2^2) q^ql + (l + 2z^) ql) + 
2M2q| (z (-3 + + (3 - 9z2) ^2^3 - z (5 + z^) q^ql + (-3 - 3^^ + 4^4^ ^3^^ ^ 

zg6 {2Ml + g|) (gjga + 4Af^ (zgi + 93)) 
128F67rgi (-4 (-1 + z2) M2 + qf) ql (4M2 + gf) ' 

^ ^ ^ (g2) gjgj (-2z^gf g3 + (1 + ql + 2A42 (2zgi + (1 + z2) gg)) ^ 

128F67r (-1 + ^2)2^2^3 

A {q,)g% (-2M2 (2zg2 + (1 + 3^2-, ^^^^ ^ 2zg|) + qs {2z\f + 2z3g2g3 ^ _ 4^2 ^ ^4-, ^^^2 _ 2zql)) 

128^677 (-l + 02)%ig3 

^(93)5^ 



(2M2 (-z2 (-3 + ^2) ^2 ^ ^ (3 + ^2) ^^^^ ^ ^ ^2) ^2) ^ 



i?5 



128^^6^ (-1 + 02)2^2^2 

93 (- + Z^) ql-{l-z^ + 2z^) qlqs + Z {l + z^) q,ql + (l + z^) ql)) - 

32F6 (-1 + .2)2 (_4 (-1 + .2) M2 + ,2) ,32 (^2^3 ((. + .) ^1 + ("l + O '^1^3 - 2zq,) + 

8(-l + .)(1 + .)M^ (3.2g? + 9z^qU3 + (-2 + 9z^ + 2z^) qiqj + z{2 + z'^) ql) + 

2MUI {z^ (-3 + z^) ql + {2z - 8z') qUs + (4 + 5z^ (-3 + z^)) q^ql + 2z (-3 + z^ + z^) ql)) + 

zg\M^ {2Ml + 932) (g2g3 + 4M2 {zq^ + gg)) 
128i^6^gi (-4 (-1 + .2) Ml + g2) ql {^Ml + gf) ' 

^ (92) (-4^2 (gi + zqz) + 93 (2z(?? + (-1 + z'') qm - '^zqj)) 



128FH {-1 + zY qiq^ 
A (93) 9\ (2M^ (z (-3 + z') - 2 (1 + z^) g^gs - 2zql) + ((l + + 2zViq3- 



USF^n {-1 + z^y qiql 
{1 + z^) qiql ~ 2zql)) + 

A {q,)!/X i:^^^ + 4^^i(j3 + (1 + :-) %) + 'y:; {-'^^<rl + (1 - 3.:-) r/i<y3 + 2zq,4 + (l + -2) ^:|)) 

128F67r (-1+^2)2^4 

I{4 : 0,-qi,q3;0)g\ ,4 //, , 2\ 2 , / , , 2\ /i , 2\ 2\ , 

Q0776^ II 2^2^ 1 , 2^ .f2 , 2^ 3 '^2 93 {{l + z ) q, + z [-1 + z ) q.q; - [l + z ) q^) + 

32F^(-l + z2) (-4 (-1 + z2) M2 + g2) ^3 

8(-l + z)(l + z)M^ {3zql + (-1 + 10z2) qfqs + 3z (l + 2z^) qiqj + (l + 2z^) ql) + 
2Mlql (z (-3 + + (3 - 9z2) ^2^3 - z (5 + z^) g,g2 ^ _ 3^2 ^ 4^4^ ^3^^ _ 

g\M^ (2M2 + ^32) (g2g3 + 4M2 {zq, + q-j)) 
128^6^ (-4 (-1 + .2) M2 + gf) ^3 (4J\^2 + ^2) ' 

^ ^(g2)gM2M^' + gi) I A{qr)g\{2z{AP^+ql)q, + q,{ml + iql + q^) 
^ 128F^TT 128F^TTqi 

A (gs) 9\ {2zqi {Ml + ^32) + (8M2 + g2 + 3^2^) 
128F67rg3 



((5 + z^) qlqlql + SM^^ (z (-3 + Az') ql+ 



128FHqr {AMI + 9?) (4 (-1 + z^) Ml - ql) 93 (4M2 + ^32) 

2 (19 - 18z2) gig3 + z (-3 + 4z2) ql) + 2M^ (4z (-1 + z^) + (77 _ 36z2) ^3^3 + 2z (33 + 8z2) ^2^2^ 
(77 - 36z2) gigf + 4z (-1 + z^) g4) + 2M2gig3 ((lO + z^) g^ + 2z (9 + 2.^) qlq^ + (29 - 7.^) g2^2^ 

2z (9 + 2z2) q^ql + (10 + z^) ql)) - 

2Ml {Aq.qs {qj + g^) + z {qf + Gqlql + qt))) , 
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^ ^ 3.gi Af^ (2Af^ + ql) 3A (93) 9\ (2M^ + gj) ((l + z^) gi + 2^93) 

^ 256F%(?2(_4(_i^^2)Af2 + g2) 256F67r(-l + z2)2g3 

3A(gi)gi(2M^ + g|) (2zgi + (l + z2)g3) ^ iAjq^) g\{2Ml + {2zql + {l + iz^) g^q^ +2zql) ^ 
2bQFH (-1 + ^2)2 ^2^3 256F67r (-1 + z-^ f qfq^ 

37(4 : 0, -gi, gg; 0)gi (2M2 + g|) 2 ^ a ^ 2n 2 ^ ^ 2n 2n 2^ ^ 
o— ;^ — ^ {—Q9 ({1 + z ) qt + z (3 + z ) O1O3 + 1 + 2; go + 

4 (-1 + z^) Ml ((1 + 2^2) (^2 + 20 (2 + q^q^ + (i + 202-) ^2^^ ^ 

^ ^zg\M^ (2M^ + gj) 3zA (gg) ff^ (2M2 + g|) ( (l + 0^) gi + 20g3) 

* 256F67rgi(-4(-l + z2)M2+g2)g3 + 256Ffi7r(-l + z2)2g2g3 

30A (gi) g6 (2M2 + g2) (2zgi + (1 + ^2) gg) (g,) g6 (2M2 + g2) (2zql + (1 + 3z2) gig3 + 20g2) 

256^6^^ (-1 + ^2)2 ^^^2 256^6,^ (-1 + Z2)2 g2g2 

3/(4:0,-gi,g3;0)zgi(2M2 + g2) 2/^,^ 2^ 2^ ^ 2n ^ ^ 2x 2x ^ 

.,^6, , , 2,2 , , , 2, ,,2 2N -«2 1 + ^ gl + ^ 3 + glg3 + 1 + g3) + 

64F6 (-1+2:2) (4 (-1 + 22)M2 - g|) g3 

4 (-1 + Ml ((1 + 2z2) g2 + 20 (2 + gigg + (1 + 2^2) ql)) , 
3^ (g2) (2M2 + g22) ((1 + 02) ql + ^ (3 + ,2) g^g3 + (1 + ,2-) ^2) 

-Kg = Ti r 

256F67r (-1 + ^2)2 ^2^2 

3^ (r/i) if'X ((1 + y") qj + 2c (2 + - [-7 + qUB + ^-tfl + ^-^-^^ ((1 + 'ji + ^^'Z^)) ^ 

256F67r(-l + z2)2g^g2 

3A (g3) g\ (20gf - 02 (_7 + ^2) g2^^ + 2z (2 + z2) g^gg ^ ^ ^2) g3 ^ 2M2 (20gi + (l + z^) gs)) ^ 

256F67r(-l + 02)^g2g3 
3/(4 : 0, -gi, g3; 0)zg\ {2Ml + gj) 2 ^ o 2 ^ / r ^ 2^ <, 2^ ^ 

64Fe(-l + 02fg,(-4(-l + 02)M2 + g2)g3 ^'^^ ^'^'^ + ^ i'^ + ^ ) ^^^S - 2g3) + 

4 (-1 + 0^) M.^ ((2 + 0^) ql + 60g.g3 + (2 + 0^) g^)) - 25aFe.g?gtf ^^^^ gi) ^ 

^ 3 (-1 + 0^) ff^M^ (2M2 + g22) 3A (g2) g% {2Ml + gj) (0gi + g3) (gi + 0g3) 
256^6^ (-4 (-1 + 02) Ml + g|) 256F67r (-1 + 02) gig3 

3A (gi) ff^ (0g3 + (1 + 202) g2g3 - (-4 + 02) g^g2 + g3 + 2M2 (0gi + g3)) 

256F67r(-l + 02)g3 

3A (g3) g\ {ql - (-4 + 02) g2g3 + (1 + 202) gigf + zgf + 2M2 (gi + 0g3)) 

2bQFH{-l + z'^)qx 
3J(4:0,-gi,g3;0)gi(2M2 + g22) / o ^ 2n ^ 2n ^ 

64Fn-l + .^)(4(-l + 02)M2-g22) (-'^^ ('^^ - - (-3 + O «1«3 + «3) + 

4 (-1 + 02) M2 ((1 + 02) g2 + 40gig3 + (l + 0^) 532)) , 



)93 



^ ^(g2)gigi(4M2 + g2+g2 
" 256^6^ (-1 + 02) g2g2 



A (g3) g^ (2M2 ((-1 + 02) g2 + 20gig3 + 2g2) + g3 (0g3 + (-1 + 202) g2g3 + zgigg + g|)) 

256F67r(_l + ^2)g2g2 

A (gi) g\ {2Ml (2g2 + 20gig3 + (-1 + 0^) g|) + gi (g? + 0g2g3 + (-1 + 202) gig| + 0gf)) 

256F67r(_l + ;j2)g2g2 

(64Fe(-1^2)^;;MM+^ll^ + g2)g2) (- (2^^ + ^) i^^l + ^3^) i^^l + + ^3^) ^ 

20^'gig3 (-4M^ + g2g32) + 02 (4M2 + g2 + ql) (4M4 + 3g2g32 + 2M2 (g2 + 532)) + 
0gig3 (8M^ + qt + qi + 4MI {ql + ql))) - 
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256i^6^g2 (4^2 + ^2) (_4 ^ ^2) ^2 + g|) (4^2 + ^f) 

AMI (2Mj + g2)%| + {2Ml + g?)^ + 2MI {AM^ + + qj) {qUI + 2M^ {qf + qj)) - 

zqiq3 (32M^ + UM^ {qf + q^) + qfqj (ql + ql) + 2M^ (qt + Aqfql + ql) ) ) . (A.2) 

In the above expressions, qi and qs are always to be understood as the magnitudes of the corresponding three-momenta 
(except in the arguments of the function /), gi = \qi \, q^ = |g3 |. Further, the function I{d : Pi,P2,P3',P4) refers to 
the scalar loop integral 



1 f 1 1 1 1 

(rf : Pi,P2,P3;Pa) - . J ^2^^rf ^ ^^^2 _M^+ie{l+ p^f -M^+ie{l+ paf - + iev {I + P4) + ie ' ^ > 



In a general case, this function depends on the four-momenta pi. For the case p^ = which we are interested in, it 
can be expressed in terms of the three-point function in Euclidean space J {d : Pi,P2,P3) 

J (d : Pi,P2,P3) = / , . — =; =; . (A.4) 

J {2nY^l+^^)2 + M^{l+P2)^+M^{l+P3)^ + M^ 
In particular, the function 7 (4 : 0, —qi, qs; 0) which enters the expressions for Ri can be written as 

/ (4 : 0, -gi, ga; 0) = i J (3 : 0, -91,9*3) . (A.5) 



For diagram (5), we obtain the following representation: 

King = Ti ■ T2 Si + ai ■ (fi(T3 • 91 Ti • T2 S'2 + CTi • 93(73 • 9i Ti • T2 6*3 + CTi • 91CT3 • 9*3 Ti ■ T2 ^4 

+ <?! • 93'?3 • 9*3 Tl ■ T2 S's + CTl • as Ti • T2 + 9l ' ^3 X <Ti Ti • T2 X T3 S7 , (A. 6) 

where the functions Si = Si{qi, 93, z) are given by 

A{q,)g\{2M^ + ql) A {q2) g\ {AM^ + qf + zqm + qI) Ajq^) g\ {2M^ + qj) 
^ USF^TT USF^TT 128FH 
J (4 : 0, -91, 93; 0) g\ {2Ml + g^) {2Ml + 9!) g\M^ 
32^6 64^6,^ ' 

g ^ A jgij tfA ((1 + <n + -^-93) A (93) tfA^i:i ('^ -'ii + (1 + '-') 1 

^ 128F67r (-1 + 22)^91 128F67r(-l + z2)2^2 

J (4 : 0, -91, 93; 0) ffl93 (-4z (-1 + z^) + 2^9f + (l + 3^^) gigs + 2^93) , 

32f6(_l + ;22)2^^ 

A (92) fli ((1 + z^) g^ + z (3 + z^) giga + (l + z^) gj) 
128F67r(-l + z2)2g2 

^ A(g3).gi((l + z^)gi+2z93) ^ ^ (gi) .gj (2zgi + (l + z^) 93) ^ (g2) (-29? + z (-5 + z^) giga - 2g|) 

^ 128F67r(-l + z2)^gi 128Ffi7r (-1 + z2)^ gs 128Ffi7r (-1 + z^)^ gig3 

J (4 : 0, -91, 93; 0) g\ (-4 (-1 + z') + (l + z^) gf + z (3 + z^) gig3 + (l + z^) gg) 

32F6(-l + z2)2 

^ ^ zA{q^)g\{{l + z^')q^+2zq^) ^ zAjq^) g\{2zqr + {l + z^) q^) ^ 
^ 128F67r(-l-Fz2)^93 12B>FH {-I + z'^f qi 

J (4 : 0, -91, 93; 0) zg\ (4z (-1 + z^) - 2zg^ - (l + 3z^) gigs - 2zgf ) 

32F6(-l + z2)^ 
zA (92) ((1 + z^) g^ + z (3 + z^) gigs + (l + z^) gj) 
128F67r(-l + z2)2 9i93 
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A(qi)g\q,{{l + z^)qi + 2zq3) A (qs) g\ {2zq, + {l + z^) qs) ^ 
USF^^n (-1 + z^f ql USF^TT (-1 + z^f q^ 

/ (4 : 0, -gi, g3; 0) gjgi (-4^ (-1 + + 2^gf + (l + 3^') gig3 + 2zg|) ^ 

32F6 (-1 + 22)2(^3 

v4 (g2) gj ((1 + ;^') gi + (3 + z') gig3 + (1 + ;^') gj) 
128i^67r(_i + z2)%| 

^(g3)glg3 (Zgl+ga) ^ (92) gj (gl - Z (-3 + g^gg + g2-) ^ ^(gpg^gl (gl +^g3) 

128i^67r(_i + 22) + 128F67r(-l + z2) 128^^67^ (_i + ^2) 

/ (4 : 0, -gl, ga; 0) g^^gigs {zqi + ga) (gi + ^^gs) 
32F6(-l + z2) 

A (gl) g\ (2M2 + 532) A {g2) gl (zgi (zgi + gg) + 2M2 (gi + zq^)) zA (gg) g\ {2M^ + g|) 



256F^TT (-1 + z2) g2 256F^TT (-1 + z2) gig2 

/ (4 : 0, -gl, ga; 0) g\ (zgi + ga) (2M2 + g2) 
64F6(-l + z2)ga 



256F67r(-l + z2)giga 



(A.7) 



Examining the above results one observes that the individual terms in the expressions for Ri and Si are singular 
for z = ±1, gl = and/or ga = 0. These singularities, however, cancel in such a way that the resulting terms in 
Eqs. (jA.ip and (jA.6P are finite. In principle, it is possible to obtain a representation for functions Ri and Si which is 
free of at least some of the singularities. In particular, the singularities at 2; = ±1 can be avoided if one expresses the 
results in terms of the functions Ji and J2 defined as 



Ji(d,gi,g3) 



'/2(rf,gl,g3) = 



1 



r{j(d:0,-gi,g-'3)-^(l + z) 



J(d:0,g*i) , J(rf:0,g-a) J(d : 0, gl + gl) 



93 + 9193 



91 + 91 93 



91 93 



2(1-) 
1 

(1 -Z2)i 



J(d:0,gl) , J(d:0,gl) , J(d:0,gl+g3) 



95 - 91 93 



91 - 91 93 
1 



91 93 



+ (1 + ^) 



:{j(d:0,-gl,gl)-^(l-z)2 



8Af2 - 2g2 + (d - l)(gi - g3)(2gi - g3) 



J(rf:0,gl) , J(d:0,gl) , J(d:0,gl+gl) 



9| - 91 93 



{d- l)g3(gi -gs)^ 



9r - 91 93 
J(rf:0,gl) 



91 93 



(A.8) 
(A.9) 



8A/2 _ 2g32 + {d- l)(gi - g3)(gi - 2g3) 



{d- l)gi(gi -g3)^ 



Jid : 0, gl) + 



2(4Af2 + (d-2)(gi-g3)^ 



(d- l)glg3(gl - g3)^ 



-Jid : 0,gi + g3) 



J(d:0,gl) , J(d:0,g*3) J(d :0,gl+gl) 



(1-z) 
8M2 - 



93+9193 9i+9i93 gig3 
8Af2 - 2g2 + {d~ l)(gi + g3)(2gi + g3) 



(d - l)g3(gl + g3)^ 
2g2 + (d-l)(gi+g3)(gi+2g3) 
(d-l)gl(gl+g3)3 



0,gl) 

2(4Af2 + (d-2)(gi+g3)^) ^ , 

^ °' (rf-l)gig3(gi+g3)^ '-'^''^ ''^\ \ /' 



rather than the three-point function J{d : 0, — gi, qa) and uses certain linear combinations of two-point functions and 
tadpole integrals. In the above expressions, the two-point function is defined as 



f d'^l 



1 



1 



(/+Pl)2+M2 (1+^-2)2+^^2 
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